Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add Burstedde Weighted BFBT test #6123

Merged
merged 1 commit into from
Nov 11, 2024

Conversation

quangx
Copy link
Contributor

@quangx quangx commented Nov 3, 2024

This pull request adds a test using weighted BFBT (with Stokes solver type = block AMG).

tests/burstedde_bfbt.sh Outdated Show resolved Hide resolved
@@ -0,0 +1,105 @@
# The 3D polynomial Stokes (burstedde) benchmark, for which an
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

can you change the comment to read "test burstedde benchmark using BFBT preconditioner", remove all the other comments here, include the burstedde.prm from the benchmarks folder and only include the relevant sections here (mesh refinement and solver settings).

@tjhei
Copy link
Member

tjhei commented Nov 8, 2024

you will need to update test results once again

@quangx quangx force-pushed the burstedde_bfbt_test branch from 764a66a to 4565bbf Compare November 9, 2024 03:11
Copy link
Member

@gassmoeller gassmoeller left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I am ok with merging this to have a test for the bfbt schur complement preconditioner, but just to document this for the future: At the moment the bfbt version requires 159 outer GMRES iterations and the simple inverse viscosity mass matrix requires only 18 outer GMRES iterations. To compute this I used the same parameter file you have here inside the tester docker container, but with set Use weighted BFBT for Schur complement = false. In other words either the method doesnt work very well for this benchmark, or there is something wrong in the method implementation. Did I interpret this correctly?

@gassmoeller gassmoeller merged commit 2bfb836 into geodynamics:main Nov 11, 2024
8 checks passed
@quangx
Copy link
Contributor Author

quangx commented Nov 12, 2024

I am ok with merging this to have a test for the bfbt schur complement preconditioner, but just to document this for the future: At the moment the bfbt version requires 159 outer GMRES iterations and the simple inverse viscosity mass matrix requires only 18 outer GMRES iterations. To compute this I used the same parameter file you have here inside the tester docker container, but with set Use weighted BFBT for Schur complement = false. In other words either the method doesnt work very well for this benchmark, or there is something wrong in the method implementation. Did I interpret this correctly?

@gassmoeller

From all the tests I have run, AMG-BFBT can take quite a lot of iterations to converge when the number of initial global refinements is small.

If you start with four initial global refinements, AMG takes 18 iterations and AMG-BFBT takes 73. For five levels,
AMG takes 19 and AMG-BFBT takes 65.

@gassmoeller
Copy link
Member

Thanks for confirming. However, is this a state you would expect? From my understanding BFBT was supposed to be a better approximation to the real Schur complement than the inverse mass matrix and therefore should reduce the number of outer GMRES iterations. Did you find a theoretical reason why it is worse?

@tjhei
Copy link
Member

tjhei commented Nov 12, 2024

Yes, we are aware. I asked him to open PRs to add tests so we can document the current status and can correctly track progress.

In short: we do not implement the corrections mentioned in the original paper for coarse meshes and I suspect our inner solver is not done correctly either.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants